On Dissipation Mechanism in the Intrinsic Josephson Effect in 
Layered Superconductors with d-wave Pairing 

S.N. Artemenko 

Institute for Radioengineering and Electronics of Russian Academy of Sciences, Moscow 103907, 

Russia 
(February 1, 2008) 

Abstract 

Conductivity mechanism in the regime of the intrinsic Josephson effect in 

layered superconductors with singlet d-wave pairing is studied theoretically. 

The cases of coherent and incoherent interlayer tunneling of electrons are 

considered. The theory with coherent tunneling describes main qualitative 

features of the effect observed in highly anisotropic high-T c superconductors 

at voltages smaller than the amplitude of the superconducting gap, mechanism 

of the resistivity being related to excitations of the quasiparticles via the d- 

wave gap. At voltages of the order of the gap amplitude and larger, the I-V 

curves are better described by the incoherent tunneling. Interaction of the 

Josephson junctions formed by the superconducting layers due to charging 

effects is shown to be small. 

PACS numbers: 74.50.+r, 74.25.Fy, 74.80.Dm 
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I. INTRODUCTION 



Theoretical understanding of the intrinsic Josephson effect (IJE) in layerd high-T c su- 
perconductors is limited by insufficient knowledge of transport properties of both the super- 
conducting and normal states. Studies of the energy structure show that electronic spectral 
density in directions (0, 7r) is affected by strong interection with spin fluctuations (see e. g. 
Ref. |T]-0). In the same time at low energies, in directions (tc, tt) corresponding to zeros of the 
order parameter A (</>), electronic structure can be described in terms of the Fermi liquid. 
This gives chances to describe the IJE at low temperatures and small voltages in terms of 
an approach based on the Fermi liquid theory. The aim of this work is to study the resistive 
properties of layered high-T c superconductors and to understand to which extent they can 
be described in the frame of the BCS model with intralayer singlet d-wave pairing. 

The IJE is expected to be qualitatively different in cases of coherent and incoherent inter- 
layer transfer of electrons. If the interlayer tunneling is coherent, a layered superconductor 
is a strongly anisotropic 3D crystal the normal-state resistivity of which is induced by scat- 
tering. On the other hand, a superconductor with incoherent tunneling can be considered 
as a stack of Josephson tunnel junctions. 

If the tunneling is incoherent the in-layer momentum is not conserved in the interlayer 
transitions. This results in a finite tunneling resistivity along the c-axis perpendicular to 
conducting layers. If the tunnel junctions are formed by conventional s-wave superconduc- 
tors the product of critical current and normal state resistivity V c = I c Rn is of the order of 
the gap value which corresponds to experimental data for high-T c superconductors. How- 
ever, for the d-wave symmetry of the order parameter, independent averaging over directions 
of the electron momentum in neighboring layers, in the case of incoherent tunneling would 
result in zero critical current along c-axis. Then in order to explain experimentally observed 
large values of critical current along c-axis one must assume a special d-wave-like angle de- 
pendence for the probability of random interlayer hopping. Fraction of the s-component in 
the order parameter in BSCCO was found in the recent tunneling measurements^ to be below 
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10~ 3 . Without assumption of the special momentum dependence of the interlayer transfer 
integral for the incoherent interlayer tunneling such symmetry would result in negligibly 
small critical current I c and in V c about three orders of magnitude smaller, than experimen- 
tally observed values^. Furthermore, with such small values of I c the regime of branching in 
which some junctions are in superconducting state while others are in a resistive state would 
be impossible in the range of voltages V ~ Ao/e per one resistive junction. Such branching 
is the most typical manifestation of the IJE. Thus it is difficult to understand main features 
of the IJE in the frame of incoherent interlayer tunneling and d-wave pairing. 

In the case of a layered superconductor with coherent tunneling, which is an anisotropic 
3D crystal, one must not assume a special angle dependence for the interlayer transfer 
integral in order to explain large experimental values of I c . But the question whether 
a clean quasi two-dimensional crystal can exhibit the IJE is not still clear. The finite 
normal-state resistivity in crystals is induced by scattering, therefore, p c oc 1/r. In a clean 
superconductor the scattering rate is small, h/r <C A . Typical voltages in experimental 
studies of IJE are by few times smaller, than the maximum energy gap A /ei, i. e., much 
larger than h/er. At such voltages (and frequencies of Josephson oscillations) resistivity 
is expected to decrease as voltage and frequency increase because of the Drude-like regime 
of scattering!, however such decaying I-V curves are not observed. The region of finite 
resistivity may become more pronounced if the intralayer scattering by impurities is resonant. 
Then gapless states are formed by such scattering at relatively large energy region around 
Fermi energy, 7 ~ Aoh/er, which results in finite conductivity at V < 7!. Calculated 
value of conductivity at voltages V <C 7 was found to be consistent with the experimental 
data! But in a clean superconductor the region where such a mechanism is expected to 
work is limited by voltages much smaller, than A . Thus, typical voltage per one junction 
in the resistive state, eV ~ A , can be explained only provided an additional mechanism 
for the resistivity is present at frequencies uj > ?i/r, 7. Such a mechanism was considered 
recently in the study of bandwidth of the Josephson plasma resonance at low temperatures!, 
it is related to the dissipation induced by electron excitation via the d-wave gap. However, 



this mechanism dies out at voltages larger than the amplitude of the d-wave gap, and to 
get finite resisitivity at higher voltages one must take into account some other mechanisms 
of scattering. A finite resistivity at large voltages, may be induced by some contribution 
of incoherent interlayer tunneling or by inelastic scattering, e.g. scattering on spin wave 
excitations is expected to become important when large energies and "hot spot" regions are 
involved. We do not consider these processes here limiting our study by low voltages. 

In this paper we study, first, the IJE in a perfect (i. e. without scattering) highly 
anisotropic superconducting crystal with coherent interlayer charge transfer and d-wave 
symmetry of the order parameter. We derive equations for charge and current densities and 
calculate I-V curves at voltages and frequencies of the order of the gap, paying attention to 
interaction between oscillations in different layers due to charging effects, considered first 
by Koyama and Tachikii. Then for comparison we calculate I-V curves for the case of 
incoherent interlayer tunneling and found that the calculated curve looks rather similar to 
the experimental data at large voltages. 



II. MAIN EQUATIONS 

Having in mind to study large frequencies and voltages much larger, than the inverse 
scattering time we calculate expressions for current and charge densities using collisionless 
transport equations in Keldysh approach for Green's functions at coinciding times G (t, t) 
derived by Volkov and Kogan0. We generalize these equations to the case of layered super- 
conductors assuming the interlayer interaction described by the tight-binding approximation 
and the superconductivity described by a BCS type Hamiltonian leading to singlet d-wave 
pairing. This approach to interlayer transitions describes the case of layered single crystal, it 
is opposite to the case of random interlayer hopping considered in Ref. [11]. Our results were 
also checked using quasiclassical transport equations for Green's functions integrated over 
momentaP generalized for layered superconductors with coherent interlayer tunneling00. 
For the case when time dependence of the scalar potential can be neglected the results can be 
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calculated using a standard expression for the conductivity in terms of the spectral densities 
in neighboring layers. Using this approach we calculate I-V curves for the case of incoherent 
interlayer tunneling. 

We consider homogeneous interlayer currents along the c-axis, this case is realized in 
narrow samples with the width in the afe-plane smaller, than the Josephson length. We solve 
equations for diagonal and off-diagonal with respect to spin indices components g n m{t) and 
fnmit) of Keldysh propagator. 
d 

o 7 9n m ^nf nm fnm^m (/^n f^m)9nm t±_ ^ ] (^-^-nn+iQn+im Qnm+i-^-rn+im) i 
at i=±l 

d 

^^777 fnm ^^fnm ^n9nm 9nm^m f/^n f^m) fnm 

at 

t± {Ann+ifn+im fnm+i-^m+im)- (X) 

i=±l 

Functions g nm (£, 4>, t), / n m(6, 4>i are matrices in layer indices. Here £ p = e (p) — ep, where 
e (p) is the single electron energy in the normal state, ep is the Fermi energy, and p is the 
electron momentum in the ab plane, is the angle of the in-plane electron momentum, t± 
is the interlayer transfer integral. Furthermore, A n = A(0) n and Xn are the anisotropic 
superconducting order parameter and its phase in layer n, /i n = e$ n + f-Jr is the gauge- 
invariant scalar potential, $ n is the electric potential, A nn+ i = exp(i(p n ), where (p n = 
Xn+i — Xn ~ i§^A z is the gauge-invariant phase difference between the layers, and A z is the 
component of the vector potential in the direction perpendicular to the layers. Electric field 
between the layers is expressed via the gauge-invariant potential as 

Ti dip n 

eE n s = fi n - Hn+i + 2 ~jp ( 2 ) 

The scalar potential ft is related to branch imbalance^l, it is responsible for charging effects 
in Josephson plasma oscillations^ and in 

The current density between layers n and n + 1 and charge density in layer n can be 
calculated as 

jn.n+1 = "Ti / /r> 4-\9 (A-n+l n9n n+1 ~ 9n+l n^-n n+l) j (3) 

2s J {2itTi) z 



Pn- U J (2^(9™ + 9* nn ), ( 4 ) 



where s is the crystalline period in c-direction. 



III. CHARGE AND CURRENT DENSITIES 

We solve equations (|l|) perturbatively with respect to t± and consider the case of low tem- 
peratures T -C A . Contributions of inter layer transitions to charge density are quadratic 
in t±, therefore, in the main approximation we can take into account only intralayer contri- 
bution which in the Fourier representation reads 

X 8A 2 /i. 



Qnn ~r 9r< 



tanh (5) 



£ e(4e 2 - uj 2 

with e = \f£ z + A 2 . Inserting this expression to (£§) we get in the time representation 



Pn 



2 

~ [ dhF^n^t-h), (6) 





where k is the inverse Thomas-Fermi screening radius, 

F(t) = /J rfe ^!|iM tanh£/2r \ (7) 

00 ' 

(• ■ ■) means averaging over <fi. 

Integral with function F describes non-exponential relaxation of fi with relaxation time 
of the order of h/A . We will need equation in the case of slow variations of the 
scalar potential when F(t) — > 25(t) so that the integral in Eq.(|5]) reduces to 2/i(t). Then 
inserting the expression for the charge density to the Poisson's equation with electric field 
determined by Eq.(|2|) we express the difference of scalar potentials between neighboring 
layers, <5/i n = fi n+ i — fi n , in terms of time derivative of phase differences. 

( / ■y 1 — j \ l m 
7S== — , (8) 
Vl + a + 1/ 

with a = 4ej_/(/€s) 2 , where e± is a high frequency dielectric constant in perpendicular direc- 
tion. For s — 15 A, 1/ n = 2 A, e± = 12 we get a ~ 0.8. Then the factor before the sum in 
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Eq.(^|) is about 0.04 and the last factor describing contributions from non-nearest neighbors 
is about 0.15. Thus, 5fi n is small compared to time derivatives of the phase differences, 
therefore, the charging effects and contribution of 5fi n to the electric field between the layers 
must be small as well (cf. Eq.(fJ)). 

To calculate the current density we need to solve equations ([I]) in the linear approximation 
in t±. But equations are still difficult to solve for arbitrary fi n (t) and <p n {t), therefore, we 
calculate expressions for current density in two limiting cases. First, we find the solutions 
in the linear approximation with respect to potential fi which describes the charging effects. 
The conditions for smallness of such effects will be discussed later. 



9n,r 



At I A 2 



tanh 



a 



2T { e(Ae 2 - ^ 



(2o;i - u)Q 



LO — UJl 



iujS, 



UJ — UJ\ 



As 2 -{to- wO 2 As 2 - lu 2 



As 2 - oj\ 



(9) 



In this limit the current density between layers n and n + 1 consists of a component deter- 
mined by the phase difference if n between the same layers only, and of a component which, 
in addition, depends on the difference of scalar potentials in these layers, 

j n>n+l {t)=f{t)+f{t). 

Assuming that t± does not depend on momentum we get 

oo 

f{t) =3cJ dtiF[ti) cos sin ; (10) 



oo OO 



f(t) = Jc J dt x J dt 2 F(t 2 ) 
o o 

<Pn{t -h- t 2 ) 



(Pnit-h) . <f n (t) 

sin ; 

2 2 

ip n (t - t x ) 



cos ■ 



cos • 



-5/i„(t - ti) 



cos • 



2 

<Pn(t) 



5fi n (t -t x - 1 2 )- 



+ sin 



8 H n (t - tx -t 2 ) sin 



<fn(t) 



11) 



2 2 2 J 

Since according to Eq.(H) 8 'f/, n (t) depends on phase differences between different layers, the 
component fllTD of the current describes interaction between the " Josephson junctions" due 
to charging effects. Interaction due to charging effects was studied in the phenomenological 
model by Koyama and Tachikii, however, the results are not similar. 

Equations (|l~0|) and ( 111) can be simplified in limiting cases. At T %u, eV Aq 



f{t) = 3c sin Vn + jc^~^{l -cos <p n ); (12) 

f{t) = 2j c jdh sin ^ ~ h) 5fi n (t - fr) sin (13) 
o 

So the current can be considered as consisting of superconducting, normal and interference 
components, and of the quasiparticle contribution related to the difference of scalar potential. 



In the limit of linear response the dissipative term in Eq.(12) vanishes because in spatially 
uniform case excitations of the quasiparticles via the superconducting gap are forbidden. In 
the nonlinear regime this dissipation mechanism is effective because in the presence of current 
along the c-axis the phase is dependent on the layer index. This makes the system non- 
uniform and excitations via the gap become allowed!. In the regime of the linear response this 
dissipation mechanism becomes possible due to scattering. Taking into account scattering 
as it was done in Ref . || we get 

2 dip n 



j(t)/jc = ¥n + 
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3rA§ dt 



2A 



The expression for current density simplifies also at large frequencies 



u, V > uj % 



(14) 



(15) 



where u p is the Josephson plasma frequency. In the most anisotropic materials like Bi- 
and Tl-based cuprates this frequency is small enough, A ^> to p . Under condition (^) the 
electronic AC current is shunted by the displacement current, 



Vac ~ V D c 



vdc, 



and time dependences of the phase differences become simple, ip n ~ 1oj n t + ip U) ^ < 1. 
When all layers are in the resistive state then 5fi = 0, and I-V curve has the form 

( V 



J 



(2A - V) 



2Ao/^ 
V V2A 



E 



9 (V - 2 A. 



K 



/2A 



E 



/2A 



V2A 
tanh 



+ 



V_ 

4T' 



(16) 



V V J V V 

where, again, V is the voltage per one junction. This expression is also valid in the limit 
a — > in which 5fi = 0. 
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At low voltages the I-V curve (|16D is described by the linear quasiparticle conductivity 
a q = nesj c / A = e 2 / (87rA 2 A ). This value agrees with the experimental data by Latyshev et 

and differs from the linear conductivity calculated for the case of the resonant scattering^ 
by factor 8/tt 2 « 1. At larger voltages the shape of the I-V curve fll^D is different from that 
observed experimentally. The calculated curve has a logarithmic peak at V = 2A , and 
current decreases with voltage at V > 2Ao, while the experimental curve si exhibit a peak in 
the conductivity at V = 2Aq and does not decay at V > 2Aq. This demonstrates that our 
approach is not valid at large V when contributions of electrons with energy of order Ao are 
important. Such electrons are strongly scattered by spin fluctuations which we did not take 
into account. One may expect that any additional mechanism of scattering, in particular, 
inelastic processes or presence of the component with the incoherent charge transfer in the 
interlayer tunneling will smear the logarithmic anomaly at V = 2Aq and add an Ohmic 
contribution to I-V curves at large voltages. 

Now we consider the regime of branching. At current values I < I c the regime is possible 
in which some "junctions" are in the superconducting state while others are in the resistive 
state. In this case the DC current through the superconducting junctions is transported in 
the form of the superconducting current, and the AC current is transported mainly as the 
displacement current. The total voltage is formed by the sum of the voltages across resistive 
junctions, and the I-V curves consist of branches corresponding to different numbers of the 
junctions in the resistive state. The number of branches is equal to the total number of 
the layers in the sample. In the limit of small a the n-th branch is described by the first 
term in Eq.(|T6D with V substituted by V/n. Under condition (|T5|) the current density was 
calculated also for arbitrary relation between 5fi n and <p n . At (p n 3> T 



]n, n+1 



0n + 
{(p n - SfXn) 3 



(17) 



2A j \ 2A 

Note that the voltage across the resistive junction differs from up n by 5fi n (cf. Eq.(^])). In the 
limit a — > 0, when 8fi n = 0, the current as function of the voltage V per one resistive junction 
is identical for all branches. At finite a, according to @, the shape of branches depends on 
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the neighboring junctions whether they are in the resistive state or not, the difference being 
smaller at smaller values of the parameter a. The shape of branches presented as function 
of the mean voltage per one junction in the resistive state V, which is defined as the total 
voltage divided by the number of resistive junctions, is shown in Fig.l and 2. One can see 
that already at a = 0.2 the branches almost coincide as it is observed experiment allySB. 
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FIGURES 




FIG. 1. I-V curve branches normalized per one layer for the system of 10 layers. The curves 
from left to right correspond to 10, 9 junctions in the resistive state, 5 and 2 neighboring junctions 
in the resistive state, 5 junctions in the resistive state separated by junctions in the superconducting 
state, and 1 junction in the resistive state. Parameter a = 0.8. Current is measured in units of J c , 
and voltage - in units of 2Aq. 




FIG. 2. Similar branches for a = 0.2. 
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IV. INCOHERENT TUNNELING 



Since a finite resistivity may be caused also by some contribution of incoherent interlayer 
tunneling at large voltages, we shall consider now briefly the case of incoherent tunneling. 

Note that expressions for the interlayer current can be derived via 2D Green's functions. 
These Green's functions are easy to calculate for the case when time dependence of /i can 
be neglected. To calculate the interlayer current we need the spectral density determined 
by the imaginary part of the Green's function. It has the form 

A(e,0 = - g*) = n[ulJ{e - ^ n + A 2 ) + vlJ(e + y/t n + A*)], (18) 
where £ n = £ + fi n , tt| n and u| n are Bogolyubov amplitudes 

u j =1(1+ ; k ] ,v 2 £ =Ml- , ^ 1 . (19) 

' n 2 v ' n 2 v v^^/ 

Then we can find the interlayer current for both coherent and incoherent tunneling. For 
the case of coherent tunneling at T = the parallel component of the momentum conserves 
in the interlayer transitions, and we find 

17/2 

j oc J dej d 2 P A(e + U/2, i n )A{e - U/2, £ n+1 ), (20) 

-U/2 

where U = (p n . We calculate this integral in the linear approximation with respect to fi 
when the result can be found easily 

i « / « ( %^ A2 H2^ - V) ) . (fifc^) . , 21 , 

This expression coinsides with the version of Eq.(|l7D linearized with respect to //. 

For the case of incoherent tunneling the momenta of the electrons in different layers are 
integrated independently, and we have 

U/2 

J oc J dej d 2 Pl A(e + U/2, Q J d 2 p 2 A(e - U/U n+1 ). (22) 

-17/2 

In this case the current is determined by the density of states in the layers 
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K ( ^ ) 9{e - A ) + — K ( — ) 6(A - e 



(23) 



V e J v ~ u/ ' A VA y 
Since the density of states does not depend on fi (at least when the latter does not depend 
on time, otherwise equation ( P^ ) is not valid), the current does not depend on Sfi as well. 
Thus the interaction between the "junctions" due to charging effects here is absent. Inserting 
Eq.(p3"D into ( p2|) we get the I-V curve presented in Fig. 3. This curve resembles experimental 
curve at large voltages, however, it does not contain the region of linear resistivity at small 
voltages. 




FIG. 3. I-V curve for the case of incoherent tunneling and all layers in the resistive state. 
Current is measured in arbitrary units, and voltage - in units of 2Aq. 



V. CONCLUSION 

Thus the model with BCS-type d-wave pairing and coherent interlayer transport de- 
scribes qualitatively such features of the IJE like branching with typical voltages per one 
junction V < Aq at low temperatures. The charging effects are found to be small at rea- 
sonable values of parameters. Damping in the system and, hence, the resistivity at voltages 
larger than inverse scattering time can be attributed to the excitations of the quasiparticles 
via the d-wave gap. At small voltages such transitions take place near the nodes of the gap, 
and the region of the angles where such transitions are possible increases with the voltage 
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increasing. The conductivity value induced by this mechanism is sample independent, in 
contrast to the damping due to scattering. 

However, our model based on the Fermi liquid approach does not describe the experi- 
mental curves at voltages of the order A and larger. Namely, in the regime when all layers 
are in the resistive state it results in the logarithmic singularity at V = A and decaying I-V 
curve at larger voltages. This discrepancy appears when contribution of the quasiparticles 
with larger energies to the conductivity becomes important, it may be related either to the 
inapplicability of the simple BSC-type model with d-wave pairing or to some other mech- 
anisms of scattering becoming effective at larger energies, like inelastic scattering on spin 
wave excitations, or other effects which smear the spectral density. We do not address these 
processes here. The I-V curve calculated for the case of incoherent tunneling looks rather 
similar to the experimental data at large voltages but is incosistent with the experimental 
data at small voltages. 

Note that at elevate temperatures, especially near T c , one may expect a different regime, 
in which the resistivity is related to quasiparticle scattering. Such mechanism is expected 
to die out at large frequencies and voltages u, V 3> 1/r. 
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